AIAA  JOURNAL,  VOL.  35,  NO.  7:  TECHNICAL  NOTES 


1235 


subsonic  turbulent  boundary  layer  is  typically  0.02 1 .  This  produced 
an  f,  =  5.34  kHz,  which  is  within  53.8%  of  Rossiter’s  predicted 
edgetones  for  the  m  =  1  mode. 

As  previously  noted,  an  examination  of  the  present  frequency 
spectra  at  0-deg  yaw  clearly  indicated22.75  kHz  as  the  fundamental 
mode,  with  secondary  peaks  between  120  and  350  Hz.  Because  the 
fundamental  mode  matches  the  acoustic  mode  f D  within  2%,  it  may 
be  suggested  that  the  dominant  oscillation  was  the  result  of  a  trans¬ 
verse  mode.  Therefore,  the  transverse  acoustic  mode  based  on  the 
cavity  depth  would  appear  to  be  the  mechanism  producingthe  pres¬ 
sure  oscillations  recorded  along  the  cavity  floor  and  not  the  mode 
correspondingto  vortex  shedding.  Yet,  Rossiter’s5  second  mode  also 
agrees  with  the  present  measurements  to  within  1.7%,  suggesting  a 
fluid  dynamic  mechanism.  In  addition,  the  second  mode  is  in  good 
agreement  with  the  computations  of  Tam  et  al., 10  who  determined 
this  frequency  numerically  via  an  TFT  of  the  pressure  history  at 
the  wall  centerline.  Their  computed  frequency  was  approximately 
26  kHz.  Further,  a  recent  work  by  the  Ref.  1 2  authors  provided  flow- 
field  simulations  that  containeda  pressure  wave  that  reflected  off  the 
floor  of  the  cavity.  This  pressure  wave  may  represent  the  transverse 
acoustic  mode. 

Although  a  combined  fluid  resonant  flow  may  be  the  logical  con¬ 
clusion,  it  is  suggested  that  the  acoustic  mode,  which  produced  a 
self-sustaining  oscillation  in  spite  of  small  geometric  changes,  may 
be  the  dominant  driving  mechanism  in  the  present  study,  at  least  un¬ 
til  yj  of  35  deg.  This  result  is  similar  to  the  flow  transition  reported 
by  Zhang  and  Edwards13  in  which  the  flow  switched  modes  from 
a  transverse  oscillation  to  a  longitudinal  oscillation  as  the  cavity 
length  was  increased.  In  the  present  study,  mode  switching  between 
3  5  and  45  deg  was  observed  to  produce  oscillations  with  frequencies 
of  either  23  or  3.5  kHz.  Beyond  45  deg,  a  permanent  mode  switch 
was  observed.  The  resultant  dominant  frequency  was  7.7  kHz  at 
45  deg,  which  increased  linearly  to  12.5  kHz  at  65  deg  yaw. 

Summary 

An  investigation  into  the  effect  of  yaw  angle  variations  in  a  two- 
dimensional  open  cavity  placed  within  a  supersonic  freestream  was 
performed.  The  effect  of  yaw  angle  changes  on  the  frequency  of 
the  dominant  pressure  oscillation  was  noted.  Specifically,  there  ap¬ 
peared  to  be  little  change  in  the  frequency  of  the  oscillation  up  to 
37.5  deg.  Comparison  of  the  measured  dominant  mode  to  the  fun¬ 
damental  acoustic  mode  based  on  cavity  depth,  i.e.,  the  transverse 
mode,  were  found  to  be  in  excellent  agreement.  This  suggests  that 
the  dominant  mechanism  for  the  present  configuration  was  acous¬ 
tic  in  nature  and  not  fluid  dynamic.  However,  a  variation  in  the 
secondary  mode,  probably  fluid  dynamic  in  nature,  showed  an  in¬ 
creased  importance  as  the  yaw  angle  was  increased. 
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Introduction 

NSTEADY  separating  flows  over  oscillating  airfoils  occur  in 
many  important  engineering  applications,  such  as  airplanes, 
helicopter  rotors  in  forward  flight,  and  turbine  blades.  There  are 
some  distinct  features  of  the  unsteady  separating  flows  (dynamic 
stall)  over  a  rapidly  oscillating  airfoil  that  draw  the  special  atten¬ 
tion  of  many  scientists.  These  features  include  large  amounts  of 
force  and  moment  hysteresis  and  oscillatory  pressure  fluctuations. 
In  most  situations,  these  features  of  dynamic  stall  significantly  limit 
the  performance  of  the  device. 

The  primary  objective  of  the  present  study  is  to  identify  the  most 
accurate,  robust,  and  economic  turbulence  model  for  dynamic  stall 
computations.  However,  testing  all  of  the  turbulence  models  avail¬ 
able  is  simply  impractical.  Alternatively,  only  a  few  popular  and  rep¬ 
resentative  models  are  selected  after  surveying  papers  and  reports 
on  separated  flows.  The  Baldwin-Lomax  (B-L)  model1  is  selected 
because  of  its  popularity  as  a  zero-equation  model.  The  Spalart- 
Allmaras  (S-A)  model2  is  chosen  among  one-equation  models  be¬ 
cause  of  its  excellent  performance.  Finally,  the  fc-e  model3  is  se¬ 
lected  because  it  is  the  most  popular  two-equation  model. 

Numerical  Methods 

The  govemingequationsare  the  time-averaged,two-dimensional, 
compressible  thin-layer  Navier-Stokes  equations.  These  governing 
equations  are  solved  by  the  TURNS  (Transonic  Unsteady  Rotor 
Navier-Stokes)  code.4  The  inviscid  convective  fluxes  E  and  G  are 
differenced  by  the  Roe’s  upwind-biased  flux-difference  splitting 
scheme.5 

The  time -marching  integrationof  the  governing  equations  is  done 
by  using  the  lower-upper  symmetric  Gauss-Seidel  (LU-SGS)  im¬ 
plicit  scheme.6  Whereas  the  LU-SGS  scheme  is  unconditionally 
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stable  and  fast,  it  is  limited  to  first-orderaccuracy  in  time.  Therefore, 
Newton  subiterations  at  each  time  step  are  used  for  misteady  flow 
computations  in  order  to  minimize  errors  arising  from  linearization 
and  factorization. 

The  turbulent  transport  equations  are  decoupled  from  the  mean 
flow  equations  and  is  solved  by  an  implicit  factored  alternating  di¬ 
rection  implicit  scheme.7 

In  dynamic  stall  computations,  a  two-dimensional, body-fitted  C- 
type  computationalgrid  moves  in  a  sinusoidalpitching  motion  about 
an  airfoil’s  quarter  chord  in  the  inertial  reference  frame  (x,  z,  t).  At 
each  time  level  the  new  location  of  the  grid  is  computed,  and  all 
metrics  and  the  Jacobian  of  the  transformation  are  recalculated.  To 
capture  the  freestreamaccurately,  a  modified  finite  volume  method8 
is  used  for  calculating  both  space  and  time  metrics. 

Results 

Unsteady  two-dimensional  turbulent  flows  over  an  oscillating 
NACA  0015  airfoil  have  been  computed  using  the  selected  three 
turbulencemodels.  A  wind-timnelexperimentconductedby  Piziali9 
at  the  U.S.  Army  Aeroflightdynamics  Directorate  has  been  chosen 
for  the  test  case. 

The  computations  were  done  for  a  freestream  Mach  number  M 
of  0.29  and  a  chord  Reynolds  number  of  1 .95  x  106.  With  a  reducea 
frequency  k  equal  to  0.1 ,  the  instantaneous an0e  of  attack  (degrees) 
of  the  pitching  airfoil  at  time  t  is  given  by 

a{t)  =  a,„  +  4.2  sin  (ikM^  +  \tc)  (1) 

where  a,,,  is  the  mean  angle  of  pitch  and  t  is  the  time  step  nondi- 
mensionalizedby  chord  and  freestream  speed  of  sound.  The  study 
considered  three  mean  angles  of  pitch  of  os„  =  4,  11,  and  15  deg, 
corresponding  to  attached  flow,  light  dynamic  stall,  and  deep  dy¬ 
namic  stall,  respectively.  The  phase  shift  \n  makes  unsteady  com¬ 
putations  begin  atthe  lowest  angle  of  attack  of  a  pitching  cycle  with 
zero  angular  velocity. 

After  investigating  the  sensitivity  of  the  computations  to  grid 
density,  time  step,  and  size  and  number  of  subiterations  at  each  time 
step,  final  computations  were  made  on  a  C-type  259  x  60  grid:  100 
points  on  the  upper  surface  and  80  points  on  the  lower  surface  of  the 
airfoil,  40  points  along  the  wake  cut,  minimum  wall  normal  spacing 
of  A:  =  0.000 1  (0.0004  for  the  k-e  model)  chord,  and  the  far-field 
boundary  located  at  1 5  chord  lengths  away  from  the  airfoil.  Time 
marching  is  achieved  using  a  time  step  size  of  At  =  0.01 13,  which 
corresponds  to  10,000  time  steps  per  pitching  cycle.  At  each  time 
step,  three  subiterations  were  used  to  reduce  numerical  errors  from 
approximate  factorization  of  the  implicit  scheme. 

Figures  1  and  2  show  the  hysteresis  of  lift,  drag,  and  pitching  mo¬ 
ment  coefficient  for  attached  flow  and  light  dynamic  stall,  respec¬ 
tively.  The  measured  force  coefficients  were  obtained  from  pressure 
data  only,  i.e.,  the  contribution  from  viscous  wall  shear  stress  was 
not  included.  In  about  half  of  a  cycle  after  the  initiation  of  unsteady 
computations,  the  hysteresis  of  force  coefficients  start  to  show  fully 
periodic  behavior  in  both  cases.  In  general,  the  computed  force  co¬ 
efficients  agree  well  with  the  experimental  data  during  the  upstroke 
for  both  flows. 

Figure  1  shows  predictions  compared  with  measurements  for  the 
attached  flow  case.  Although  the  flow  is  simple,  there  is  enough 
boundary-  layer  thickening  to  make  it  difficult  to  predict  the  hys¬ 
teresis  of  the  force  coefficients  accurately.  The  predictions  by  the 
S-A  and  the  k-e  models  agree  very  well  with  the  measured  data 
for  all  three  force  coefficients.  The  B  L  model  shows  fairly  good 
agreement  with  the  measurements  for  C;  and  Cj  but  not  for  Cm . 

In  light  dynamic  stall,  shown  in  Fig.  2,  the  flow  separates  near  the 
trailing  edge  as  the  airfoil  pitches  up,  and  this  separation  disappears 
as  the  airfoil  pitches  down.  Because  the  light  dynamic  stall  is  sen¬ 
sitive  to  the  onset  of  separation,  it  is  a  good  test  case  for  evaluating 
turbulence  models.  The  B-L  model  performs  poorly  in  predicting 
all  three  force  coefficients.  The  thin  envelope  of  the  C/  hysteresis 
curve  by  the  B-L  model  indicates  that  the  model  does  not  give  sep¬ 
aration  during  the  oscillatory  cycle.  The  k-e  model  shows  excellent 
agreement  with  the  measured  C/  whereas  the  S-A  model  shows 
a  little  underprediction  during  the  downstroke.  For  the  drag  and 
moment  coefficients,  the  S-A  and  the  k-e  models  show  comparable 
agreement  with  the  experimental  data. 


Fig.  1  Hysteresis  of  lift,  drag,  and  pitching  moment  for  attached  flow 
of  an  oscillating  NACA  0015  airfoil  at  (X  (deg)  =  4  +  4.2  sin(2t M  t  + 

1.5tt),  Mm  =  0.29,  Re  =  1.95  x  106,  k  =  0.1: - ,  S-A; - ,  B-L; 

- ,  k-e ;  and  a  experiment. 


angle 


Fig.  2  Hysteresis  of  lift,  drag,  and  pitching  moment  for  light  dynamic 
stall  of  an  oscillating  NACA  0015  airfoil  at  (X  (deg)  =  11+4.2  sin(2A  JVf  f 

+  1.57T),  Mm  =  0.29,  Re  =  1.95  x  106,  k  =  0.1: - ,  S-A; - ,  B-L; 

- ,  k-e ;  and  a  experiment. 
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Fig.  4  Vorticity  contours  for  oscillating  NACA  0015  airfoil:  (X  (deg) 
=  15  +  4.2sin(2tMoof  +  1.5tt),  Mm  =  0.29,  Re  =  1.95  x  106,  and 


The  computed  deep-stallhysteresis  loops  of  the  force  coefficients 
(not  shown)  are  less  satisfactory,  and  they  can  only  be  considered 
semiquantitativeon  the  downstroke.  However,  Fig.  3  shows  a  series 
of  instantaneous  pictures  of  computed  streaklines  of  the  deep  dy¬ 
namic  stall  case.  During  the  second  cycle  of  computation  with  the 
S-A  model,  solutions  and  grids  at  600  intermediate  time  steps  were 
saved  for  computing  the  streaklines.  The  computation  was  done  by 
the  unsteady  flow  analysis  toolkit  (UFAT). 10  At  each  time  step,  new 
particles  are  released  from  two  near-leading-edge  stations:  one  on 
the  upper  and  the  other  on  the  lower  surface  of  the  airfoil.  Then 
UFAT  follows  these  particles,  revealing  a  Lagrangian  description  of 
flow  phenomena. 

Figure  4  shows  a  series  of  calculated  vorticity  contours  for  the 
deepdynamic  stall  case.  According  to  the  calculations, the  vorticity 
of  the  flowfield  varies  widely  fromzero  to  a  few  thousand. The  large- 
scale  vorticesobserved  in  Fig.  3  have  low  vorticity  in  the  order  of  one 
while  vorticity  greater  than  1 0  is  limited  to  the  wall  region  near  the 
leading  edge.  Therefore,  the  vorticity  contour  levels  ranging  from  0 
to  10  is  used.  In  this  way,  the  evolution  of  large-scale  vortices  can 
be  captured  without  losing  the  global  picture  of  the  flow  structure. 
The  vorticity  near  the  leading  edge  of  the  upper  surface  is  the  most 
intense  due  to  the  flow  acceleration  along  the  curved  leading  edge 
of  the  upper  surface.  During  the  whole  cycle,  the  vorticity  contours 
give  the  overall  appearance  of  remaining  attached  to  the  leading 
and  trailing  edges.  Only  their  size  and  thicknesses  vary  as  the  airfoil 
pitches  up  and  down. 


Fig.3  Particle  tracing  of  oscillating  NACA0015  airfoil:  (X  (deg)  =  15  + 
4.2sin(24Moot+  1.5tt),  =  0.29,  Re=  1.95  x  106,  and  k  =  0.1. 
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Nomenclature 

e  -  specific  energy,  Cv  T 
f  =  flow  parameters  (p,  U,  V,  T) 

P  -  pressure 

0  =  heat  flux 

Re  =  Reynolds  number 

T  =  temperature 

t  =  time 


ters  governs  heat  flux  distribution  along  the  surface  (along  the  X 
coordinate): 

(PocUo6To^  QAX 2),  QJXf]  (1) 

We  seek  to  determine  p^.  U ^  and  T  using  heat  flux  measure¬ 
ments  Qw(  Xj)  as  the  input  data.  We  shall  minimize  the  discrepancy 
between  computed  and  experimental  values  of  Q,,.(  Xj),  i.e., 

^00  UOC  Tot  =  2“p(^)  -  Qc:mp(Kj)f  (2) 

The  mathematical  model  consists  of  a  flowfield  solver  and  an 
optimization  algorithm  minimizing  the  discrepancy  of  computed 
and  experimentaldata  by  varying  the  boundary  conditions. 

The  Jacoby  matrix  A,j  =  3g(X,)/3/}  for  the  transformation 
(p,  T,  U ) >Q(X 1,  X2,  X3)  is  studied  for  the  singularity  search: 

d£K£ 1)  dQjx 1)  dQ(X{)~ 

3  p  3  U  3  T 

dQjXA  3  Q(X2)  3  Q(X2) 

dp  3 U  dT  (  ’ 

d_Q{X 3)  3  Q(  £3)  dQ(X2) 

dp  dU  dT 

The  determinant  is  identically  equal  to  zero  in  the  following  events 
(global  singularity).  1)  The  flowfield  is  determined  not  by  the  total 
parameter  space  (p_  t/,  T)  but  by  some  subspace  of  the  lower 
dimension,  for  example,  (Re^  2)  The  heat  flux  Qw(t,  X)  is 
determined  not  by  total  space  of  flow  parameters  but  by  some  sub¬ 
space  of  lowerdimension.  3)  The  spatialheat  flux  variation  QK(X)  is 
multiplicative  along  the  X coordinate:  QW(X)=  Qn,(X0)F(X/  X0), 
i.e.,  the  heat  flux  is  self-similar  along  X. 

The  computation  of  the  matrix  A  determinant  allows  one  to 
consider  the  singularity  only  locally  in  regard  to  the  flow  pa¬ 
rameters  because  of  the  problem  nonlinearity.  A  global  singular¬ 
ity  can  be  detected  by  means  of  the  Lie  group.  In  particular,  the 
transformation  Q(X)  =  Q(X0,  p ^  T^tflXl  X0)  corresponds 

to  the  Navier-Stokes  equations  invariance  relative  the  operator 
B  =  TfPoQ  Uqq  7.^3/  3  X.  This  operator  is  absent  in  the  set  of  the 
Lie  group  basis  oithe  Navier-Stokes  equations  that  provides  the 


SQiX,) 

~  Vj  ~ 


U  —  velocity 

X,  Y  =  coordinates 

p  =  viscosity 


absence  of  this  type  of  singularity.  Nevertheless,  the  transformation 
Q{X)  =  Q{X0,  p^  TqApIX/  Xo)  can  take  the  place  for  an 
infinitely  thin  boundary  MyenThus,  although  the  matrix  (3)  is  non¬ 
singular,  it  can  be  very  close  to  a  singular  one  that  causes  instability. 


Subscripts  and  Superscripts 

comp  =  computation 
exp  =  experiment 
w  =  wall  conditions 

=  freestream  parameters 

Introduction 

HEN  it  is  necessary  to  determine  flow  parameters  govern¬ 
ing  heat  transfer,  the  need  to  solve  inverse  convection  prob¬ 
lems  arises.  For  example,  determining  the  boundary  heat  flux  us¬ 
ing  temperature  measurements  at  another  boundary  is  discussed  by 
Moutsoglou1  for  a  two-dimensional  region  with  steady-state  flow 
describedby  the  Navier-Stokes  equations.  Estimation  of  freestream 
parameters  is  of  practical  importance  in  many  problems .  Direct  mea¬ 
surement  of  these  parameters  is  accompanied  by  certain  technical 
difficulties  caused,  for  example,  by  flowfield  distortion  or  the  heat 
flux  acting  on  the  sensor.  Wall  temperature  measurements  are  the 
simplest  from  this  viewpoint. 

The  estimation  of  free  flow  parameters  using  heat  flux  measure¬ 
ments  at  different  points  at  body  surface  is  considered.  Let  us  con¬ 
sider  the  supersonic  gas  flow  described  by  equations  of  viscous 
compressible  gas  (Navier-Stokes).  The  set  of  freestream  parame¬ 


A  global  singularity  exists  for  special  events  and  can  be  easily 
detected  by  computations.  The  most  practical  interesting  problem 
arises  when  a  global  singularity  is  absent.  In  that  event,  the  non¬ 
linearity  of  the  problem  becomes  crucial.  Qw(t,  X)  nonlinearly  de¬ 
pends  on  (Pqq  U^,  T)  and  nonlinearly  varies  along  X.  Then,  there 

is  a  six  parameterlocal  transformation:  (X3  y/3)  P _ >Q3-  If 

this  transformation  is  smooth  and  if  p  andX>3  are  smooth  mani¬ 
folds,  then,  according  to  the  Sard  theorem,2  such  a  transformation 
should  be  nonsingular  (in  the  common  case).  If  det[dQ( X,)/ 3/;] 
is  not  identical  to  zero,  the  set  of  points  X  p  R3  for  which 
det[3  Q(  Xj )/  dfj  ]  =  0  is  either  empty  or  a  smoothtwo-dimensional 
surface  embedded  into  R3  space.  It  can  be  easily  seen  that  f3  and  Q3 
are  smooth  manifolds,  and  the  transformation  is  smooth  for  flows 
considered  here  (having  no  shocks  on  the  surface).  Thus,  the  prob¬ 
lem  [Eq.  (1)]  will  be  nonsingular  with  the  probability  close  to  1. 

These  speculations  lead  to  the  conclusion  that,  if  the  numerical 
experiments  have  demonstrated  nonsingularity  in  some  point,  the 
problem  is  nonsingular  in  a  significant  part  of  the  phase  space. 

A  number  of  difficulties  arise  at  the  investigation  of  the  Navier- 
Stokes  problem,  in  general,  which  are  connected  with  the  need  for 
repeated  solution  of  the  direct  problem  and  large  computation  time. 
In  this  connection,  the  two-dimensional  spatially  evolutional  flows 
that  can  be  solved  with  quick  marching  methods  are  of  special  inter¬ 
est.  The  supersonic  flow  over  a  flat  plate  with  a  compression  corner 
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